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ADIO  INTERFERENCE  DUE  TO  ROCKET  EXHAUST  JETS: 

THE  MEASUREMENT  AND  COMPUTER  MODELLING  OF  AMPLITUDE  AND 

PHASE  NOISE  SPECTRAj 


by 


Williams 


SUMMARY 


A  full  description  is  given  of  a  computer  model  developed  to  describe 
amplitude  and  phase  noise  spectra  generated  by  incoherent  volume  scattering  in 
the  turbulent  ionised  exhaust  jet. 

The  experimental  method  and  equipment  are  described  in  detail  and,  by 
choice  of  a  computer  model  for  the  jet  velocity  and  density  structure,  a 
comparison  is  made  between  the  measured  and  the  computed  noise  power  spectrum 
density  functions.  The  agreement  is  excellent  for  the  particular  moderately 
ionised  jet  selected.  Velocity  turbulence  is  not  found  to  be  significant  in 
establishing  the  noise  spectrum.  The  radial  and  longitudinal  gradients  in  the 
steady  mean  local  velocity  within  the  jet  adequately  predict  the  observed 
spectrum  function. 


Paper  presented  at  the  Sixth  JANNAF  Conference  on  Rocket  Exhaust  Plume 
Technology,  Monterey,  California,  March  1971 
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1  INTRODUCTION 

The  addition  of  amplitude  noise  to  a  microwave  carrier  by  turbulent, 
ionised  rocket  exhaust  jets  has  been  studied  in  the  UK  and  spectra  recorded 
since  the  early  nineteen  sixties,  A  system  capable  of  measuring  both  phase 
and  amplitude  noise  powers  was  added  in  1967  and  sited  at  the  RPE,  Westcott. 

The  system  permits  noise  spectra  to  be  recorded  down  to  -120  db/Hz  below  the 
carrier  level  and  out  to  200  KHz  off  the  carrier  frequency. 

In  1966  a  simple  model  was  proposed^  for  radio  wave  scattering  in  rocket 

exhaust  jets  and  found  capable  of  explaining  much  of  the  gross  character  of 

2 

the  amplitude  noise  spectra  observed.  Early  computations  agreed  well  with 
experiment . 

It  is  the  purpose  of  this  report  to  describe  in  some  detail  the  form  of 
the  computer  model  now  in  use,  to  give  a  brief  description  of  the  microwave 
instrumentation,  and  finally  to  compare  computed  noise  power  spectrum  functions 
with  the  experimental  data. 

In  constructing  a  quantitative  model  of  electromagnetic  wave  scattering  by 
a  turbulent  rocket  exhaust  jet  it  is  necessary  to  make  physical  and 
mathematical  approximations  and  assumptions.  This  is  because  of  both  the 
complexity  of  the  problem  and  the  lack  of  knowledge  in  a  number  of  important 
areas . 

To  test  the  scattering  model  a  detailed  knowledge  of  jet  structure  is 
required.  Furthermore,  accurate  and  unambiguous  experimental  data  are  only 
available  for  rocket  motors  tethered  at  ground  level.  As  a  result  the  jet 
description  chosen  is  particular  to  such  motors.  The  essential  property  of 

this  model  is  simplicity  resulting  in  ease  of  computation.  There  are  other 

3  4b  ...  3  . 

models  available  5  9  of  more  general  applicability,  one  of  which  is 

regularly  in  use  at  the  RPE,  and  there  are  certain  to  be  other  and  better  ones 
in  the  future.  The  overall  method  of  approach  has  therefore  been  one  of 
flexibility  throughout.  In  particular  the  algebraic  union  of  jet  and  electro¬ 
magnetic  model  has  been  avoided.  Quadrature  is  used  to  evaluate  a  volume 
integral  for  the  received  signal  due  to  scattering  over  the  whole  jet.  Thus 
any  concept  of  an  exhaust  jet  can  be  considered,  the  jet  being  input  to  the 
scattering  program  as  a  number  of  two-dimensional  ( ax i symmetric )  arrays  from 
a  preceding  computer  program. 

2  ELEMENTS  OF  THE  MODEL 

A  diagrammatic  description  of  radio  wave  single  scattering  by  a  local 
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turbulent  volume  in  a  rocket  motor  exhaust  jet  is  shown  in  Fig.  1.  A  small 
fraction  of  the  total  scattered  energy  enters  the  receiving  antenna  at  a 
frequency  Doppler  shifted  with  respect  to  the  incident  frequency  by  virtue  of 
the  relative  motion  of  the  scattering  element  dV  with  respect  to  the  antennae„ 
When  all  such  volumes  within  the  jet  are  summed  the  result  is  a  spread  of  energy 
about  the  direct  received  carrier  wave  frequency.  The  magnitude  and  shape  of 
this  noise  spectrum  is  the  desired  end  product. 

Neglecting  losses  along  the  incident  and  scattered  ray  paths,  but 
remembering  that  they  must  be  included  in  the  final  calculation,  we  may  write 
the  scattered  power  density  at  the  receiver  in  the  form 


PT  CT  (<V 

U  7T  r  2 


a  dV 


U  7T 


where  a  is  the  volume  scattering  cross  section  per  unit  volume  defined  such 
that  if  the  total  power  contained  in  a  section  of  the  incident  wave  front  having 
an  area  a  dV  were  radiated  by  an  isotropic  radiator  located  at  the  element  dV, 
the  strength  of  the  wave  reaching  the  receiving  antenna  would  be  the  same  as 
the  actual  strength  produced  by  scattering  in  the  volume  element  dV  (for 
other  symbols  see  Nomenclature) . 

The  received  scattered  power  is 


dP, 


ft  gt  jV 

It  It  r  2 


7T  r 


x  s  (v 

U  7T 


dV 


The  total  noise  power  received  is  therefore 


(h  TT  ) 


3  T 


; 


gt  (ai}  gr  ((V 

vol  r  j'  r22 


a  dV 


where  the  volume  of  integration  includes  the  jet. 
The  direct  received  signal  power  is 


P  =  P 
S  T 


Gt  (a2)  Gr  (B2) 


(H  tt) 
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In  the  limit  of  a  transmitter  effectively  at  infinity 


^  "  Gr  (62) 


gr  (61> 


o  dV 


vol 


(1) 


The  Doppler  frequency  shift  of  the  power  scattered  from  the  element  dV  is 


f 


D 


u  / 

—  (cos 


cos  B  ) 


or,  for  the  transmitter  at  infinity, 


fD  =  (cos  e2  -  cos  ) 


(2) 


This  description  has  for  simplicity  been  two  dimensional  onlyc  It  is 
readily  extended  to  three  dimensions  for  computation  purposes.  In  fact  for 
all  the  experimental  data  available  the  transmitting  antenna  lies  in  the  plane 
defined  by  the  receiving  antenna  and  the  jet  axis. 

3  THE  VOLUME  SCATTERING  CROSS  SECTION 

The  derivation  of  the  differential  scattering  cross  section  for  a  plane 
wave  incident  on  a  turbulent,  under-dense  plasma  medium  is  given  by  Tatar ski 
(Chapter  U) . 

In  the  limit  of  the  first  Born  approximation  and  following  an  electro¬ 
magnetic  analysis  he  obtains 


da  =  2  7T  k  sin  4>  (2  k  sin  ^*)  dQ  dV 


It  follows  that  the  cross  section  per  unit  volume  is 


a 


k^  sin2  \p  4>n(2  k  sin 


(3) 
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The  spectrum  function  <f>^  of  the  refractive  index  fluctuations  is  not 
known  with  any  accuracy  at  the  moment.  Salpeter  and  Treiman^  in  a  useful 
paper  discuss  some  of  its  general  properties  and  characteristics,  listing 
amongst  the  possibilities  the  spectrum  function  based  on  the  Kolmogoroff 
theory  of  turbulence  previously  discussed  in  general  terms  by  Tatarski  ,  who 
gives 


t>n  (2  k  sin  — ) 


r  (n  +  f  ) 
ir  /it  r  ( n ) 


2  3 

n  a 


,  ,  ,  2,2  .  2  y 

1  +  4  a  k  sm 


n  +  3/, 


Putting  n  =  —  , 


((>  (2  k  sin  ?-)  = 

n  2 


0.062  n  2  a3 


,  ,  2,2  .  2  v 

I  +  4  a  k  sm  -7 


"/, 


( **) 


To  a  good  approximation  at  microwave  frequencies 


n 


2  ti  N  eL 

2  2 

m  ( a)  +  v  ) 


With  N 
index 


N  +  N1  we  obtain  for  the  mean  square  fluctuation  in  refractive 


h 

~ 

m 


(w‘  + 


2.2 

v  ) 


(5) 


Substitution  gives 


k  2 

32  it  r 


a  = 


^77  I2  (N)2 


2  3 

sin  \\)  0.062  a 


1  + 


.2 


2  2  2  Y 

1  +  )i  a  k  sin  77 


11/, 


cm 


-1 


(6) 


where  I  = 


(n^2) 


is  the  turbulent  intensity. 


CONFIDENTIAL 


7 


CONFIDENTIAL 


The  Born  solution  (equation  (3))  assumes  that  the  incident  wave  is 
everywhere  unperturbed  by  the  turbulent  medium.  In  fact  as  the  electron 
density  increases  both  the  incident  and  scattered  wave  will  suffer  attenuation 
through  absorption  and  scattering  loss.  In  this  model  an  attempt  is  made  to 
extend  the  validity  of  the  single  scatter  description  to  attenuating  media  by 
calculation  of  these  losses.  Along  both  the  incident  and  scattered  ray  the 
absorption  per  unit  path  length  is  taken  as 


6 

a 


0.46  N  v 

2  2 

0)  +  v 


db  cm 


The  loss  due  to  scattering  along  incident  and  scattered  rays  is 


6 

s 


a 

s 


0.2303 


db  cm 


where  a  is  the  total  scattering  cross  section*  To  obtain  this  the 
s 

differential  scattering  cross  section  must  be  integrated  over  all  solid  angles. 
We  obtain 


a 

s 


1  p  p  —2  o 

8  tt  r  IN  0*062  a 
_ e _ 

2 ,  2n2 
( 1  +  v  /u)  ) 


Q 


where 


Q 


6tt 

3 

P 


-  j  (2p?  +  2p  +  1)  {  ( 1  +  2p)  5/6 


1}  - 


-  2(1  +  p)  {  (1  +  2p) 1 /6  -  1  }  +  y  {  (1  +  2p)  7/6  -  1  } 


and  p 


2  a 


O 


The  total  attenuation  per  unit  path  length  is  6  =  S  +  S 

.Las 

Integration  of  6^  along  both  incident  and  scattered  paths  to  and  from  a 
volume  element  dV  ,  taking  account  of  the  variation  in  plasma  properties 
along  the  paths,  gives  the  reduction  in  the  received  scattered  power  caused 
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by  absorption  and  volume  scatteringo  This  may  be  evaluated  readily  if  all 
plasma  properties  are  known. 

)|  THK  ROCKET  JET  MODEL10 

Several  models  of  turbulent  reacting  flow  are  already  in  existence  5  * 

A  rigorous  treatment  of  the  physical  and  chemical  processes  involved  in  the 
mixing  of  a  hot  combustible  gas  with  air  has  not  been  attempted  in  this  work; 
rather  the  emphasis  has  been  placed  on  a  mainly  empirical  approach  which  allows 
a  rapid  and  simple  estimate  of  the  gross  features  of  the  jet.  The  treatment 
falls  naturally  and  conveniently  into  two  parts;  the  fluid  flow  field  of  the 
jet  is  described  first,  followed  by  a  description  of  the  electron  density 
field . 

The  method  used  to  determine  the  nozzle  exit  plane  condition,  which 
constitutes  the  starting  point  of  the  jet  model,  is  based  on  an  implicit,  one¬ 
dimensional  approximation  of  nozzle  flow,  including  finite  rate  chemistry. 

4 . 1  The  jet  flow  field 

The  jet  consists  of  a  mixing  region  of  annular  cross  section  and  a  core 
consisting  only  of  exhaust  products;  the  annulus  increases  in  width 
approximately  linearly  with  distance  from  the  exit  plane  until  at  some  axial 
distance  L  it  extends  across  the  whole  width  of  the  jet.  The  distance  L 
is  termed  the  core  length,  because  it  defines  the  axial  extent  of  the  core 
of  undisturbed  jet  fluid,  roughly  conical  in  shape,  based  on  the  nozzle  exit 
diameter.  Beyond  the  core,  where  the  jet  gas  is  mixing  with  the  atmosphere 
at  all  points  within  the  jet,  two  regions  have  been  described  by  experimental 
observers:  first  a  transition  region,  over  which  the  effect  of  the  core  and 
any  particular  flow  properties  are  dissipated,  and  secondly  what  is  termed 
the  fully  developed  region,  in  which  it  is  possible  to  find  semi-empirical 
formulae  for  many  of  the  jet  properties.  This  region  -  the  "far  field"  of 
the  fluid  flow  -  is  usually  the  least  interesting  insofar  as  radio 
interference  is  concerned.  However  use  is  made  of  the  expressions  derived  for 
this  region  as  in  rocket  exhaust  jets  the  transition  region  is  relatively 
short.  Because  of  this  and  the  difficulty  of  finding  an  adequate  description 
of  the  transition  region  it  is  assumed  in  the  jet  model  that  transition  from 
the  core  region  to  the  fully  developed  mixing  region  takes  place  over  a 
negligibly  small  distance. 

The  conical  nozzle,  which  has  an  expanding  section  in  the  form  of  a 
truncated  cone,  is  the  most  widely  used  nozzle  design.  There  are  no  conditions 
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under  which  the  conical  nozzle  will  produce  a  jet  which  is  free  of  shocks. 

Even  when  the  exit  plane  pressure  is  equal  to  the  ambient  pressure  the  approx¬ 
imately  radial  flow  produced  by  the  expansion  cone  means  that  the  jet  gas 
continues  to  expand  after  it  has  left  the  nozzle.  A  full  description  of  the 
jet  must  take  some  account  of  the  shocks  formed  within  the  flow.  However, 

Q 

experimental  evidence  has  been  presented  for  non-reacting  jets  to  suggest 
that  to  a  first  approximation  the  shock  structure  does  not  affect  the  flow 
field  of  the  jet.  Hence  it  is  assumed  that  the  jet  may  be  equated  to  the  jet 
from  a  correctly  expanded  contoured  nozzle,  for  which  the  central  core  has 
constant  properties  throughout  its  length. 

A  correction  is  applied  to  allow  for  the  pressure  adjustment  in  the  jet 

which  occurs  when  the  motor  is  operated  at  other  than  its  design  conditions  <> 

The  jet  gas  is  allowed  to  expand  inviscidly  from  the  nozzle  lip  until  the 

tangent  to  the  inviscid  jet  boundary  is  parallel  to  the  jet  axis,  Fig.  2. 

The  calculation  of  the  inviscid  portion  of  the  jet  when  the  exit  plane 

pressure  is  greater  than  the  ambient  pressure  is  accomplished  by  means  of  a 

11 

method  of  characteristics  computer  program  of  general  applicability  ;  when 
the  exit  plane  pressure  is  less  than  that  in  the  surrounding  air  the 
curvature  of  the  shock  and  jet  boundary  are  determined  by  an  approximate 
method.  This  method  of  approximating  the  separating  streamline  is  set  out  in 
the  Appendix.  Having  established  the  initial  portion  of  the  jet  it  is  assumed 
that  mixing  takes  place  from  this  point  onwards. 

This  procedure  is  a  convenient,  if  much  simplified,  way  of  allowing  for 
the  under-  or  over-expansion.  Its  major  effect  is  one  of  changing  the  scale 
of  the  jet;  that  part  of  the  jet  downstream  of  x^  (Eig.  2)  behaves  as  if  it 
were  produced  by  a  parallel-flow  nozzle  of  radius  r^  .  The  core  and  mixing 
layer  are  accordingly  non-dimens ionalized  by  the  radius  r^  .  This  approach 
to  the  problem  of  incorrect  expansion  is  similar  to  that  suggested  by  Donaldson 
and  Gray  .  It  is  acknowledged  that  this  simple  treatment  will  not  adequately 
describe  the  mixing  region  over  the  first  portion  of  the  jet;  indeed 
comparisons  of  calculated  and  experimental  attenuation  levels  close  to  the 
jet  exit  plane  have  indicated,  as  might  be  expected,  that  some  allowance  for 
mixing  must  be  made  in  this  section  of  the  flow.  However,  the  success  of  this 
correction  in  describing  the  overall  scale  and  general  behaviour  of  the  rest 
of  the  jet  is  such  that  no  major  modification  is  envisaged  at  the  time  of 
writing. 

The  model  is  made  quantitative  by  the  introduction  of  several  empirical 
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relationships.  The  core  length  L  is  assumed  to  be  a  function  of  the  Mach 

number  M  on  the  axis  at  x  ,  and  is  given  bv 
m  m  5  D  47 


L 


r 


m 


2.1  M  2 
m 


(T) 


This  simple  relationship  provides  a  reliable  estimate  of  the  longitudinal 
scale  of  the  jet,  and  has  been  tested  for  several  practical  cases  under  widely 
different  exit  plane  conditions.  Typical  results  are  tabulated  below: 


v-/v 

.1  m 


L/r 


m 


equation  (7) 


L,  calculated, 

metres 


L,  observed, 

metres 


Rocket 

1 

0.07 

21.5 

0.23 

0.22 

Rocket 

2 

0.5 

25. ^ 

2.1 

1  .9 

Rocket 

3 

2.0 

32.5 

3.7 

U.O 

The  core  is  assumed  to  be  a  cone,  centred  on  the  jet  axis,  of  base  radius  r 

*  m 

and  height  L  ,  so  that  the  generator  is  given  by 


r . 

l 

r 

m 


=  1 


X 

L  # 


(8) 


The  gas  velocity  is  assumed  to  remain  constant  throughout  the  core,  and  is 
allowed  to  decay  on  the  axis  beyond  the  core  inversely  with  x  ,  such  that 


u  x 

core 


(9) 


The  rate  of  decay  of  velocity  on  the  jet  axis  has  been  the  subject  of 

experimental  studies  (see,  for  example,  ref.  9)  in  which  the  observed  rate  has 

-1  -2 

been  measured  variously  between  x  and  x  .  Measurements  made  at  the 
RPE  in  cold  non-reacting  jets  have  indicated  x  ^  and  this  dependence  is  used 
in  the  model  in  the  absence  of  any  demonstrably  better  relationship  . 

The  radial  variation  of  the  longitudinal  velocity  component  in  the  mixing 


CONFIDENTIAL 


CONFIDENTIAL 


1  1 


region  is  expressed  as  an  exponential  function  of  position 


—  =  exp 

u. 

1 


-  ln(2) 


(  1 
r  -  r . 

l 

ro 

_ 1 

r i  -  r. 

2  i 

J 

(10) 


where 


u. 

i 

=  u  if 

core 

X 

< 

L  ; 

u . 

1 

=  u 

if 

X 

> 

L  ; 

r . 

l 

=  0 

if 

X 

> 

L  ; 

ri 

is  the 

radial 

co- 

-ordinate 

u./2 

i 
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This  expression  is  similar  to  that  used  by  Donaldson  and  Gray  but  has 
been  modified  to  give  a  continuous  value  of  <Su/6r  at  r  =  r^  „  The  use  of 
the  same  function  for  both  the  annular  mixing  region  and  the  fully  developed 
mixing  region  implies  that  the  whole  of  the  turbulent  mixing  layer  is  self- 
preserving.  Although  self-preservation  of  the  velocity  profile  is  a  property 
of  the  fully  developed  region  it  is  not  strictly  applicable  to  the  rest  of  the 
mixing  region.  However  measurements  of  velocity  profiles  in  a  rocket  jet  have 
established  that  the  use  of  equation  (10)  gives  sufficiently  good  agreement 
with  experimental  results. 

To  establish  an  expression  for  ri  it  is  necessary  to  consider  the 

2 

longitudinal  momentum  in  the  jet*  The  integral  of  the  rate  of  flow  of 
momentum  across  a  plane  perpendicular  to  the  axis  is  a  constant  of  the  system 


p  u  2i\r  dr  =  constant  = 


2  2 

p  .  u  .  Trr  . 
J  J  J 


(id 


in  which  p.  is  the  exit  plane  density 
J 

u.  is  the  exit  plane  velocity 
J 

r.  is  the  exit  plane  radius 
J 

As  u  is  a  known  integrable  function  of  r  (equation  (10))  it  is  only 

necessary  to  define  p  as  a  function  of  r  to  be  able  to  deduce  ri  as  a 

2 

function  of  position.  In  a  cold  non-reacting  jet  it  is  possible  to  invoke 
the  similarity  relationship. 
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p  -  p 


Pi  '  Pa 


Sc. 


u 

u. 


(12) 


where  Sc  is  the  turbulent  Schmidt  number  and  subscript  a  refers  to 
ambient  conditions.  For  a  hot  combusting  rocket  exhaust  jet  where  temperatures 
are  very  much  higher  than  the  ambient  temperature  over  most  of  the  jet  it  is 
sufficient  to  assume  the  approximation 


p(r)  =  constant  =  p 

core 


03) 


The  error  introduced  into  the  momentum  integral  by  this  approximation  may  be 
neglected,  as  at  the  edge  of  the  jet,  when 


P 


p  >  p 

a  core 


the  velocity  has  fallen  to  a  small  fraction  of  u_.  and  the  contribution  to 
the  integral  in  this  part  of  the  jet  is  small. 

From  equations  (6)  and  (13),  with  the  further  assumption  that 


2  12  2 

u  r  dr  =  —  u .  r .  , 

2  J  J 

o 


r . 

f  1  2 

u  r  dr 
'  o 


But  for  0  <  r  <  r.  we  have  u  =  u  =  u.  so  that  equation  (1^)  becomes 

l  core  j 

2 

r°°  U . 

u^  r  dr  =  7^-  (r?  +  r?)  .  (15) 

Jr.  c  J  x 


2  , 

u  r  dr 


r . 
i 


1  2  2 

—  u.  r. 

2  J  J 


0*0 
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Substitution  for  u  from  equation  (10)  gives 


r°° 

f  2 

u  exp 

-  2  ln( 2) 

r  r  -  r.  ^ 
i 

“ 

ip  ^ 

"  L  J 

- 

r  dr 


2 

U1  2  2 

P*1  (ri  +  rp  .  (16) 


Writing  r  =  r  -  r^  ,  dr  =  dr  ,  the  integral  is  transformed  to 


u!  (  r  +  r . )  exp 
J  1 


o 


where 


so  that 


2 

z 


2  In  (2) 

( r !  -  r.  ) 
2  1 


/rr 


2 

u.  r . 

l  l  2z 


(17) 


(18) 


Equation  (  1 8 )  is  a  quadratic  in  r-,  .  If  x  <;  L 

2 

he  found  in  terms  of  r  ,  which  is  a  function  of 

If  x  >  L  then  r.  =  0  and  u.  =  u.  L/x  so  that 

1  i  J 

as  a  function  of  x 


then  u.  =  u.  and  ri  can 
1  J  2 

x  as  given  by  equation  (8). 

rn  can  be  found  directly 
2 


This  defines  the  complete  flow  field. 

4.2  The  electron  density  field 

A  proper  determination  of  the  electron  density  distribution  throughout 
the  jet  requires  a  knowledge  of  the  gas  temperature  and  composition  and,  if 
the  full  rigour  of  finite  rate  chemistry  is  to  be  applied,  of  the  chemical 
reaction  rates  and  of  a  detailed  history  of  each  part  of  the  flow  field.  The 
complexity  which  an  exact  solution  would  produce  is  undesirable  within  the 
terms  of  reference  of  this  exercise.  A  compromise  has  been  reached  therefore 
between  the  limits  of  frozen  and  equilibrium  chemistry,  at  least  so  far  as 
the  electron  density  distribution  is  concerned. 
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If  the  chemistry  of  the  jet,  downstream  of  the  nozzle  exit  plane,  is 
assumed  to  be  frozen  (i.e.  vanishing  reaction  rates)  there  can  be  no 
recombination,  no  secondary  combustion,  no  heat  release:  the  jet  must  behave 
in  all  respects  like  a  single  jet  of  heated  inert  gas*  Under  this  condition 
the  fraction  f  of  initial  jet  gas  per  unit  mass  of  gas  in  the  mixing  region 
is  described  by  observed  empirical  relationships.  The  decay  of  jet  gas 

concentration  on  the  jet  axis  beyond  the  end  of  the  core  has  been  found  to 

”2 

vary  as  x  (ref.  9)  and  the  radial  profile  of  concentration  is  governed  by 
the  similarity  relationship 


f 
f . 


u 

u. 


(19) 


As  in  equation  (  1 0 ) ,  the  subscript  refers  to  core  values  if  x  ^  L  and  centre¬ 
line  values  if  x  >  L  •  The  first  step  in  determining  the  electron  density 
distribution  in  the  model  is  to  assume  that  the  mixing  region  geometry  is 
defined  by  the  velocity  flow  field  described  in  the  previous  section  and  that 
the  (pre-combustion)  composition  of  the  gas  mixture  is  described  by  the 
variation  of  f  with  position.  Here  Scf  is  assumed  to  be  unity,  as  there  is 
no  better  estimate  appropriate  to  combusting  systems o  This  drastic  simplifi¬ 
cation  of  the  mixing  process  is  tantamount  to  assuming  instantaneous  mixing  of 
the  jet  gas  within  the  confines  of  the  velocity  field.  The  value  of  f  at  any 
point  is  given  by 


f 

f 


1  inside  the  core 


f\  exp 


ln(2) 


f . 
1 


x  $  L 


x  >  L 


/ 


(20) 


Having  in  this  way  specified  a  mixture  composition  at  all  points  within  the 
mixing  region  the  effects  of  secondary  combustion  and  heat  release  are 
simulated  by  allowing  instantaneous  reaction  to  equilibrium,  making  proper 
allowance  for  the  velocity  profile  in  the  energy  balance  equation.  The 
ionization  reactions  are  assumed  to  be  sufficiently  rapid  in  the  hottest  parts 
of  the  mixing  region  to  maintain  equilibrium,,  Similarly  in  the  outer  edge  of 
the  jet  the  gas  velocity  is  relatively  low,  so  that  the  residence  time  of  a 
mass  of  gas  within  the  range  of  interest  of  the  flow  field  (about  two  core 
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lengths)  is  fairly  long;  it  is  assumed  to  be  sufficiently  long  to  allow  the 
ionic  reactions  to  maintain  equilibrium.  Equilibrium  electron  densities  are 
therefore  assumed  at  the  outer  edge  of  the  jet  and  inwards  to  about  the  centre 
of  the  mixing  region.  Close  to  the  central  core  and,  beyond  the  core,  near 
the  jet  axis,  the  equilibrium  electron  levels  are  much  lower  than  would  be  the 
case  if  the  frozen  chemistry  approximation  were  applied.  It  is  known  that  the 
exit  plane  electron  density,  which  is  used  throughout  the  core,  is  well  above 
the  equilibrium  concentration.  Accordingly  in  the  inner  part  of  the  mixing 
region  the  electron  density  follows  the  profile  produced  by  turbulent  diffusion 
from  the  core,  i.e*  equation  (19).  The  resultant  electron  density  profile  is 
typically  illustrated  by  Fig.  3.  In  detail  the  electron  density  at  a  radial 
point  r  is  given  by 


unless 


when 


N  =  f (r )  N 

core 


f(r)  N 


< 

core 


N 

equil 


(r) 


N 


N  ... 
equil 


(r) 


; 


(21) 


1  o  -3 

For  the  motor  to  be  described  N  =  2„10  cm  ,  and  the  maximum  N 

10-3  C°re 

in  the  mixing  layer  is  8.10  cm 

No  apologies  are  made  for  the  lack  of  rigour  in  deriving  this  model:  its 
justifications  are  its  relative  simplicity,  ease  of  application  and,  under  sea- 
level  conditions  at  least,  its  reproduction  of  what  are  believed  to  be  the  major 
features  of  the  electron  density  profile  in  a  fuel-rich  rocket  exhaust  jet,  i.e. 
the  relatively  low  centre-line  value  with  a  highly  ionised  annular  region. 

4.3  Collision  frequency 

The  electron-neutral  collision  frequency  at  any  point  is  calculated  from 


v  =  v  Z  N.  q. 

e  .  l  i 
l 


(22) 


where 

v 

e 

N. 

l 

qi 


1 8k  T 

7T  m 


is  the  electron  velocity,  with  usual  notation; 

is  the  number  density  of  species  i  ; 

is  the  electron  collision  cross  section  of  species 


CONFIDENTIAL 


CONFIDENTIAL 


16 


In  practice  only  a  few  of  the  species  present  in  the  exhaust  gases  have 

collision  cross  sections  large  enough  to  contribute  to  the  collision  frequency „ 

These  are  typically  Ho0  ,  ,  CO  and  CO  .  Thus  the  calculation  of 

the  equilibrium  composition  and  temperature  in  the  mixing  region,  described  in 

the  preparation  of  the  electron  density  profile,  serves  also  as  the  basis  for 

determining  collision  frequency,  v  is  found  to  vary  little  throughout  the 

11  -1 

jet.  It  has  therefore  been  kept  constant  at  a  value  of  2.10  sec 
h.h  The  turbulent  scale 

The  evaluation  of  the  scattering  cross  section  per  unit  volume  requires 

a  knowledge  of  the  turbulent  scale  a  (equation  (6)).  The  particular  scale 

of  interest  is  that  appropriate  to  the  fluctuating  electron  density  field  in 

a  chemically  reacting  turbulent  flow.  It  has  been  shown  both  experimentally 
12 

and  theoretically  that  this  scale  may  be  dependent  upon  the  detailed  rates 
of  electron  removal  (recombination)  in  laboratory  turbulent  diffusion  flames. 

h 

In  particular  at  a  Reynolds  number  of  ^  10  in  a  flame  of  hydrogen/nitrogen 
combusting  in  air  it  has  been  found  that  transferring  the  dominant  electron 
recombination  process  from 

H^O  +  e  neutral  products 

to 

+  — 

N  a  +  e  +  M-^N  +  M 

a 

increases  the  scale  by  a  factor  greater  than  two.  According  to  the  spectrum 

function  4>  this  could  make  a  difference  of  -10  db  to  the  power  scattered 
n 

at  small  angles. 

Perhaps  all  that  can  be  said  at  this  time  is  that,  at  least  toward  the 

edge  of  the  jet  where  the  "local"  Reynolds  number  is  low,  finite  rate  ionic 

chemistry  could  well  play  a  part  in  determining  the  turbulent  scale.  For  our 

purpose  we  shall  ignore  this  possibility.  We  shall  assume  the  electron 

density  to  be  in  the  nature  of  a  passive  additive,  frozen  into  the  flow.  This 

.  13 

has  certainly  been  shown  to  approximate  to  the  truth  by  Granatstem  and  by 

.  1  )|  .  . 

Garosi  in  the  high  Reynolds  number  flow  of  weakly  ionised  argon.  The  latter 

finds  a  regime  in  the  turbulent  ion  density  power  spectrum  associated  with  the 

inertial  subrange  of  Kolmogoroff  in  agreement  with  our  choice  for  <f>^ 

As  the  mixing/shear  layer  grows  in  width  downstream  of  the  motor  exit, 
so  the  turbulent  scale  must  grow  in  size  as  the  turbulence  develops. 
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In  a  recently  published  work  Fisher  and  Johnston^  tentatively 
conclude  that  the  three-dimensional  power  spectrum  in  a  Mach  3.3U  cold  jet 
approximates  the  -11/3  power  law,  but  more  usefully  they  determine  the 
dependence  of  turbulent  scale  on  shear  layer  width.  For  the  angles  of 
illumination  of  interest  to  this  study  we  are  dependent  upon  the  radial  and 
tangential  scales  defined  by  these  authors  in  terms  of 


■  oo 

l  U  ,  o)  d£ 


where  i  (£  ,  t)  is  the  cross  correlation  coefficient  between  two  points  in 
the  flow  separated  by  the  distance  £  and  delayed  by  the  time  interval  x 
As  may  be  seen,  this  scale  is  twice  the  scale  normally  used. 

From  Fig.  32  of  Fisher  and  Johnston, 


L 

y 


w 

3 


or 


a 


To  a  good  approximation  in  the  jet  model  used,  W  the  width  of  the 
mixing  layer  is  given  by 

2  r  x 

w  =  - 2L_ 

L 

or 


a 


r,  x 
m 

3  L 


(23) 


Thus  at  the  end  of  the  core  (x  =  L)  the  turbulent  scale  equals  r  /3  which 

m 

for  the  motor  considered  here  is  approximately  3  cm.  On  the  assumption  that 
the  electron  density  in  our  case  and  the  concentration  of  "tracer"  in  that  of 
Fisher  and  Johnston  are  both  passive  additives  in  the  flow  and  that  there  are 
no  loss  mechanisms  in  either  which  have  effective  rates  faster  than  that  due  to 
turbulent  mixing,  the  identification  of  the  two  scales  as  equal  in  magnitude 
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is  legitimate. 

The  expression  (23)  is  therefore  used  throughout  the  jet  to  determine 
the  local  turbulent  scale. 

h .  5  Turbulent  intensity 

The  turbulent  intensity  is  difficult  to  determine  in  the  absence  of 
experimental  measurements  of  the  fluctuating  electron  density  field  in  the 
rocket  jet„  However  there  is  a  fair  body  of  evidence  characterised  by 
Granat stein  and  Buchsbaum^'  that  at  high  Reynolds  numbers  the  intensity  may 
be  taken  as  unity.  This  is  accepted  and  applied  throughout  the  jet.  The 
possible  error  introduced  by  such  an  assumption  may  be  assessed  at  perhaps 
6  db  at  the  most  when  due  account  is  taken  of  the  dependence  of  scattering 
loss  on  this  quantity  and  of  the  generally  falling  mean  electron  density  with 
decreasing  intensity  at  the  edges  of  the  jet.  This  is  not  considered  to  be 
too  significant  within  the  measuring  accuracy  of  ±  3  db. 

5  THE  METHOD  OF  COMPUTATION 

The  method  selected  is  entirely  numerical,  being  based  upon  the  division 
of  the  space  comprising  the  jet  into  finite  size  volume  elements. 

The  receiver  is  assumed  to  be  on  the  "vehicle  body",  in  the  exit  plane 
of  the  motor  and  a  distance  "l"  from  the  jet  axis.  The  jet  is  defined  on  a 
cylindrical  co-ordinate  system  (x,  r,  0)  where  x  is  the  distance  from  the 
exit  plane,  r  is  the  radius  and  0  is  the  azimuthal  angle.  The  size  of  the 
volume  element  dV  =  r  dx  dr  d9  is  determined  by  an  accuracy  factor  g  .  The 
number  of  elements  of  each  lamina  is  given  at  any  station  in  the  jet  by 


n 


r 


g 


ne  =  2  6 


n 


x 


g  X 

max 

R 

max 


where  X 

max 

lengths  and 


is  the  effective  end  of  the  jet,  invariably  taken  as  two  cone 

R  is  the  maximum  radius  of  the  jet,  defined  by  the 

max 


—  =  0.  1  contour 
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It  has  been  found  by  trial  calculations  that  for  an  accuracy  of  ±1  db  in  the 
computed  noise  spectrum  the  total  number  of  elements  must  exceed  10*\ 

As  the  computer  steps  through  these  volume  elements,  considering  every 
one  in  turn,  it  determines  the  velocity  at  the  center  of  the  element  by  inter¬ 
polating  within  a  two-dimensional  (axisymmetric )  array  of  local  jet  velocities 
previously  calculated  by  a  data  production  program  based  on  Section  h .  The 
condition  for  inclusion  of  an  individual  element  is  that 


where  f^  is  the  Doppler  frequecy  shift  and 

bw 

2 

bw 

2  * 

f  is  one  of  a  predetermined  number  of  frequency  ’bins1  of  bandwidth  bw  , 
all  required  f  and  bw  being  entered  as  input. 

If  a  particular  element  is  to  be  included,  then  the  contribution  to  the 
noise  power  in  the  frequency  bin  labelled  Tf^T  is  W  r  dr  d0  dx  ,  less  the 
attenuation,  where  W  is  the  "weighting  factor’1  of  the  element  of  volume 
r  dr  d0  dx  given  by 


W  = 


^  ^  2 

32  7T  r 
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0.062  a" 


2- 

’♦*2 

00  ■ 
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1  +  h  a  k  sm  ■*- 
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It  has  been  found  convenient  to  enter  the  experimentally  determined  receiving 
antenna  gain  as  a  polynomial  in  terms  of  cos  0/2  ,  in  accordance  with  a 
least  squares  fit  to  the  experimental  data* 

The  attenuation  for  each  element  is  dependent  upon  its  position  in  the 

jet  and  the  actual  paths  taken  by  the  incident  and  scattered  rays.  Using  the 

previously  given  total  attenuation  per  unit  path  length  (6  =6  +  6  )  the 

j.  a  s 

method  of  calculating  the  total  loss  is  to  proceed  in  certain  step  lengths 
towards  the  receiver  as  far  as  the  edge  of  the  jet  and  then  towards  the  trans¬ 
mitter  as  far  as  the  edge  of  the  jet  evaluating  the  line  integral 
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Attenuation 


ST(s)  ds 
s 


along  both  paths. 

The  step  length  is  itself  a  function  of  position  in  the  jet  given  by 


N . 
1 


N 


L 

q. 


where  L  is  the  length  of  the  potential  core  and  q  is  h  g  .  This  permits 
larger  step  lengths  to  be  considered  as  the  electron  density,  and  hence  loss 
per  unit  step  length,  decreases. 

Direct  attenuation,  in  a  straight  line  from  receiver  to  transmitter  is 
also  computed,  see  Table  1. 

20 

6  EXPERIMENTAL  EQUIPMENT 

6. 1  The  amplitude  and  phase  noise  measuring  system 

The  equipment  described  here  was  designed  and  manufactured  by 
GEC  (Electronics)  Ltd  to  an  agreed  specification  as  a  result  of  a  feasibility 
study  conducted  by  this  same  company. 

The  basic  system  consists  of  an  X-band  microwave  bridge  in  which  one  arm 
is  an  air  path  subjected  to  modulation  by  rocket  exhausts  and  the  second  arm  is 
a  coaxial  cable  link  providing  an  unmodulated  reference. 

The  transmission  path  lies  between  two  aerials  31  m  apart  and  mounted  on 
top  of  towers  9  m  high.  A  cabin  is  situated  half-way  up  each  tower,  one 
containing  the  transmitter  and  the  other  the  input  stages  of  the  receiver. 

The  electronic  equipment  in  each  tower  is  protected  by  a  specially  constructed 
enclosure  from  any  excessive  acoustic  vibration  that  may  occur. 

A  control  room  placed  90  m  from  the  towers  houses  the  Main  Assembly. 

This  equipment  processes  the  signals  from  the  receiver,  extracts  the  amplitude 
and  phase  modulation  components  and  converts  them  into  a  form  suitable  for 
recording  on  tape.  Interconnections  between  the  three  Assemblies  enable  the 
complete  system  to  be  controlled  from  the  Main  Assembly.  A  block  diagram  of 
the  system  is  shown  in  Fig.  hm 

The  transmitter  is  a  klystron  operating  at  X-band  and  supplying  an  output 
power  of  1.U  watts  into  a  ,,hoghornn  aerial.  An  additional  low  power  output 
feeds  the  reference  path  which  is  a  coaxial  link  between  the  transmitter  and 
receiver  assemblies. 
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The  receiving  aerial  is  adjacent  to  the  rocket  exhaust  and  the  microwave 
signal  is  passed  to  a  mixer  where  it  is  mixed  with  the  local  oscillator  signal 
to  produce  an  intermediate  frequency  of  20  MHz.  A  second  mixer  processes  the 
reference  signal  in  a  similar  manner.  The  two  IF  signals  are  then  amplified 
and  coupled  back  to  the  Main  Assembly  through  coaxial  links. 

Further  amplification  takes  place  before  diode  detection  is  used  to 
separate  the  amplitude  modulation  from  the  carrier.  The  IF  signals  are  also 
compared  in  a  phase  sensitive  detector  to  obtain  the  phase  modulation  levels. 
The  attenuation  produced  by  the  rocket  efflux  is  measured  by  extracting  part 
of  the  input  to  the  signal  channel  IF  amplifier  and  feeding  it  to  a  separate 
amplifier  with  manual  gain  control. 

Following  each  modulation  detector  are  three  high  pass  filters  and 
associated  amplifiers  enabling  the  overall  modulation  spectrum  to  be  divided 
into  discrete  frequency  bands  suitable  for  recording  on  tape. 

To  check  that  physical  movement  of  the  receiver  aerial  does  not  introduce 
spurious  sidebands,  three  transducers  are  mounted  on  the  aerial,  one  in  each 
plane. 


Performance 


Attenuation 

Dynamic  range  -U0  db 

System  Resolution  Amplitude  and  phase  modulation 


0  db  RF  attenuation 

30  db  RF  attenuation 

Cross  Talk 

Signal  AM  to  PM 
PM  to  Signal  AM 
Ref  AM  to  PM 


-100  db  to  -120  db 
-99  db  to  -103  db 


-27  db 
-30  db 
-10  db 


6.2  Calibration 

Overall  calibration  of  the  system  is  performed  by  extracting  a  small, 
precisely  known  amount  of  the  transmitter  output  signal,  modulating  it  100$ 
in  amplitude  with  a  square  wave  and  returning  it  to  the  main  path  prior  to 
transmission.  Both  amplitude  and  phase  modulation  may  be  obtained  in  this  way 
by  returning  the  modulated  signal  either  in  phase  or  quadrature  with  the 
carrier.  The  small  sample  of  transmitter  power  is  coupled  into  the  auxiliary 
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calibration  circuit  from  a  20  db  cross  coupler,  via  a  set  level  attenuator, 
into  a  3  db  sidewall  coupler.  This  coupler  divides  the  available  power 
equally  between  two  identical  branches  wherein  it  is  subjected  to  100# 
modulation  by  means  of  diode  modulators  and  recombined  in  a  second  3  db 
sidewall  coupler.  The  modulated  wave  is  then  reintroduced  into  the  main 
waveguide  run  via  another  20  db  cross  coupler  and  adjusted,  by  means  of  a 
phase  shifter  in  each  branch,  to  be  either  in  phase  for  amplitude  modulation 
or  quadrature  for  phase  modulation.  Ideally,  this  system  requires  that  the 
modulator  should  produce  side  bands  only,  which  requires  a  balanced  modulator 
and  careful  setting  up.  However,  where  the  sideband  is  low,  the  level  of  the 
reintroduced  carrier  is  also  low  and  is  negligible  in  comparison  with  the 
carrier  provided  by  the  unattenuated  transmitter  output.  An  effect  is  only 
apparent  in  the  case  of  phase  modulation,  where  it  produces  an  amplitude 
modulated  component,  in  this  case  more  than  Uo  db  down  on  the  phase  modulation 
calibration. 

6.3  Receiver  aerials 

Of  particular  significance  to  the  measurement  of  small  levels  of  phase 
noise  is  the  mounting  of  the  antenna  in  the  vicinity  of  the  rocket  jet,  in  this 
case  the  receiving  antenna.  The  magnitude  of  the  problem  is  illustrated  by 
the  fact  that  on  an  X-band  carrier  a  phase  modulation  of  -130  db^  is  produced 
by  a  physical  movement  of  the  antenna  by  one  microinch«  As  a  consequence, 
following  a  study  of  the  hostile  environment,  a  receiving  aerial  was  designed 
supported  by  six  rubber  shock  mounts  in  a  steel  housing.  This  is  shown 
mounted  alongside  the  motor  in  Fig.  5. 

The  feed  from  the  aerial  consists  of  a  short  length  of  coaxitube  to 
provide  a  flexible  non-phase  shifting  connection. 

Vibration  transducers  are  used  during  all  firings  to  determine  the 
motion  of  the  receiving  aerial  in  three  dimensions.  Interpretation  of  these 
data  allows  a  full  measurement  of  that  part  of  the  phase  noise  spectrum,  if 
any,  which  is  attributable  to  acoustical  vibration.  In  the  results  reported 
here  this  contribution  proved  to  be  negligible  at  all  frequencies. 

6.U  Recording  and  analysis 

Data  concentrated  in  the  Main  Assembly  are  brought  to  a  lU-channel 
Ampex  FR  1300  tape  recorder.  The  information  recorded  is  listed  in  Table  2. 

It  will  be  noted  that  channel  capacity  is  divided  between  flame  interference 
data  and  a  check  on  system  performance.  Both  direct  recording  (DR)  and 
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frequency  modulated  (FM)  modules  are  employed  giving  a  frequency  recording 
capability  from  direct  voltage  to  300  kHz  at  a  tape  speed  of  60  inches  per 
second.  The  specified  response  "flatness*’  is  to  within  ±  3  db;  this  is  easy 
to  obtain  and  can  be  improved  with  careful  adjustment  during  alignment.  A 
dynamic  range  of  between  25  db  and  32  db  is  available  on  DR  channels  depending 
upon  tape  speed,  whilst  that  for  FM  channels  is  between  40  db  and  kh  db,  again 
dependent  upon  tape  speed,  the  total  harmonic  distortion  not  exceeding  2.0%. 

Attenuation  levels  may  also  be  displayed  on  a  pen  recorder  simultaneously 
with  recording  on  tape,  thus  permitting  instantaneous  information  regarding 
this  variable  to  be  obtained. 

Analysis  is  performed  using  the  James  Scott  Spectrum  Analyser  Type  190 
GHLA  and  the  BrUel  and  Kjaer  Level  Recorder  Type  2305.  A  block  diagram  is 
shown  in  Fig.  6. 

The  spectrum  analyser  makes  a  Fourier  analysis  of  the  noise  power 
available  from  the  amplitude  and  phase  modulation  channels  of  the  tape  recorder. 
Controlled  manually,  the  levels  of  noise  input  at  any  frequency  between  100  Hz 
and  200  kHz  can  be  obtained,  measured  in  a  pre-selected  bandwidth.  Available 
bandwidths  are  30  Hz,  70  Hz,  100  Hz,  500  Hz  and  1000  Hz  (Broadband).  It  will 
be  noticed  that  two  attachments  form  part  of  the  spectrum  analyser,  a  Low 
Frequency  Attachment  (LFA)  and  a  High  Frequency  Attachment  (HFA)e  These  serve 
to  extend  the  frequency  coverage  of  the  main  unit  which  is  limited  to  frequency 
operation  in  the  range  1.0  kHz  to  1^5  kHz.  The  LFA  extends  the  lower  frequency 
range  to  100  Hz  whilst  the  HFA  increases  the  upper  range  to  200  kHz.  The 
intrinsic  noise  of  the  analyser  becomes  apparent  at  -120  db  with  respect  to 
1.0  volt  when  measured  in  a  100  Hz  bandwidth,  and  an  accuracy  of  ±  3  db  is 
achieved  throughout  its  frequency  range. 

For  analysis  presentation  a  Bruel  and  Kjaer  level  recorder  is  used. 

This  instrument  accepts  signals  in  the  frequency  range  2  Hz  to  200  kHz  when 
driven  in  the  AC  mode  (three  options:  rms ,  average  or  peak).  When  direct 
voltage  is  applied  it  is  converted  to  a  square  wave  by  an  electromechanical 
"chopper"  operating  at  100  Hz.  The  subsequent  process  is  then  identical  to 
that  for  alternating  voltages.  The  dynamic  range  of  the  recorder  is 
determined  by  interchangeable  range  potentiometers  having  logarithmic  or  linear 
characteristics . 

7  THE  ROCKET  MOTOR 

The  motor  used  during  these  trials  was  the  ABL  EM-71  motor  containing 
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i*0  lb  of  ELP  propellent.  It  is  well  described  in  ref.  17.  (Note  that  the 
EM-71  is  a  heavyweight  version  of  the  EM-27.  They  only  differ  in  chamber 
wall  thickness  and  end  plate  design.)  Suffice  to  say  here  that  the  propellent 
was  composite-modified  double-base  of  the  following  composition: 


Ingredient 


Per  cent, 
by  weight 


Nitrocellulose  ( 12*6/5  N)  28.0 

Nitroglycerin  38.3 

Ammonium  perchlorate  19.6 

Mg  (30)/Al  (70)  alloy  k .9 

Triacetin  7.2 

Resorcinol  1 o0 

2-nitrodiphenylamine  1 .0 

Alkali  metal  content  ppm 

Sodium  U6.5 

Potassium  10 

Caesium  5 


The  propellent  burned  at  a  calculated  temperature  of  about  3200°K  and 
at  a  chamber  pressure  of  260  psi,  the  combustion  products  being  expanded 
through  a  nozzle  of  expansion  ratio  5.^.  The  nozzle  exit  plane  diameter  was 
6  inches.  The  fuel  index  of  the  exhaust  gases  was  approximately  0oU,  and  the 
exit  gas  temperature  approximately  1700°K.  A  thrust  of  2500  lb  was  developed 
for  a  duration  just  less  than  *4  seconds. 

8  RESULTS  AND  DISCUSSION 

The  table  gives  a  succinct  comparison  between  measured  insertion  loss 
and  computed  attenuation  of  the  direct  ray  versus  angle  of  look  a  .  The 
experimental  data  have  been  averaged  over  the  firing  period  and  over  a  number 
of  firings,  and  have  a  spread  at  any  angle  of  at  least  ±  1  db. 

Comparison  of  computed  and  experimental  insertion  loss 


Angle  of  look, 

a 

Insertion 

Measured 

loss,  db 
Computed 

+  5° 

5 

7=6 

+9° 

6.5 

6.5 

+30° 

2 

1 .2 
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The  agreement  is  good.  This  is  taken  as  adequate  verification  of  the  jet 

1 8 

model  described  here.  As  shown  elsewhere  the  dependence  of  the  longitudinal 
insertion  loss  on  angle  of  look  is  a  searching  test  of  the  adequacy  of  any 
description  of  the  ionised  jet  in  the  absence  of  diffraction* 

^ig.  7  and  8  show  the  experimental  amplitude  and  phase  noise  power 
spectra,  the  bar  indicating  the  likely  error.  The  vertical  scale  is  in  terms 
of  db  below  the  incident  carrier  per  unit  bandwidth,  the  incident  carrier  being 
defined  as  the  received  signal  level  in  the  absence  of  the  exhaust  jet.  The 
levelling  off  of  the  spectra  at  'v  -120  db  at  the  highest  frequencies  is  due  to 
overal  system  noise  -  particularly  the  magnetic-tape  noise. 

The  first  thing  to  notice  is  that  amplitude  and  phase  noise  spectra  are 

identical  under  all  the  experimental  conditions.  This  is  a  first  requirement 

for  the  scattering  model  to  be  correct.  It  may  readily  be  shown  that  the 

incoherent  summation  of  the  scattered  power  from  each  correlated  volume  in  the 

jet  produces  identical  amplitude  and  phase  noise  spectra  when  quoted  in  db 

below  incident  carrier.  The  condition  is  that  the  received  scattered  power  is 

much  less  than  the  received  carrier  power.  The  agreement  between  computed  and 

measured  spectra  is  excellent.  Particularly  noteworthy  is  the  internal 

consistency  in  the  results o  At  large  angles  of  look  (±  30°)  the  low  frequency 

saturation  is  we] 1  reproduced  by  the  computed  data.  The  approximate  10  db 

difference  between  the  low  frequency  levels  at  +30°  and  -30°  is  well  reproduced. 

In  fact  the  computer  output  follows  the  changes  in  shape  and  in  magnitude  of 

the  measured  noise  spectra  over  the  full  range  of  the  experimental  data. 

However  the  originally  computed  spectra  are  below  those  measured  at  all  angles, 

particularly  at  the  large  ones*  The  very  nature  of  the  scattering  cross  section 

is  such  that  most  power  is  scattered  about  the  forward  direction;  scattering  at 

large  angles  is  weak.  At  small  angles  the  scattering  cross  section  varies  as 

the  cube  of  the  turbulent  scale.  One  is  therefore  tempted  to  increase  the 

scales  beyond  those  derived  from  the  work  of  Fisher  and  Johnston.  The  effect 

of  multiplying  by  a  factor  of  four  is  shown  as  the  broken  line  in  Fig.  7  and 

0.  The  agreement  between  experiment  and  model  at  positive  angles  is  improved 

in  this  way  but,  as  one  might  expect,  at  large  negative  angles  the  opposite 

is  true.  The  reason  is  simple  in  that  at  these  negative  angles  the  angle  of 

scatter  into  the  receiver  is  necessarily  large  and  at  large  angles  of  scatter 

-2/3 

the  scattered  power  varies  as  a  .  Numerical  trials  have  shown  that  the 

relationship  chosen  to  describe  the  turbulent  scale  gives  the  best  possible 
fit  to  the  experimental  data. 
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It  is  now  with  confidence  that  one  may  perform  a  parametric  study 
covering  the  many  other  variables  in  the  model.  Although  such  a  study  has 
been  completed,  a  full  discussion  here  would  extend  this  paper  unduly.  A 
separate  report  is  intended. 

Finally,  it  is  only  fair  to  comment  that  success  has  been  achieved  in 
the  endeavour  to  descibe  numerically  the  amplitude  and  phase  noise,  with  the 
simplest  of  scattering  models,  the  validity  of  which  may  be  questioned  on 
certain  counts.  For  example  no  account  has  been  taken  of  the  velocity 
turbulence.  In  principle  the  random  nature  of  the  magnitude  and  direction 
of  the  velocity  vector  must  spread  the  frequency  content  of  the  scattered 
signal.  In  this  model  the  frequency  spread  in  the  noise  spectrum  stems  solely 
from  gradients  in  the  steady  mean  velocity  flow  field.  An  elemental  volume 
dV  ,  of  linear  dimension  less  than  the  scale  of  this  gradient,  scatters  a 
monochromatic  signal  into  the  receiving  antenna.  With  the  inclusion  of 
velocity  turbulence  this  signal  must  be  spread  across  the  spectrum. 

A  further  paper  will  describe  a  Monte  Carlo  simulation  of  the 
scattering  problem.  In  this  simulation  all  events  are  described  in  terms  of 
probability  distributions.  Far  ranging  flexibility  is  attained  at  the  expense 
of  computing  time.  In  particular  one  has  the  option  of  adding  a  random 
component  to  the  mean  velocity  vector,  described  in  terms  of  a  normal  distribu¬ 
tion  with  zero  mean  and  with  standard  deviation  as  input 0  The  result  of 
including  an  isotropic  velocity  turbulence  of  standard  deviation  (intensity) 

O.U  is  shown  in  Fig.  9  for  the  case  a  =  -30°  .  Runs  were  completed  at 

a  =  +  5°  and  +30°  but  in  these  no  significant  effect  of  velocity  turbulence 
can  be  discerned. 

As  plotted  in  Fig.  9  velocity  turbulence  apparently  lifts  the  power 
spectrum  at  the  highest  frequencies.  (Not  shown  is  the  fact  that  this  also 
happens  at  the  lowest  frequencies  (<  100  Hz).)  At  first  sight  an  improved 
agreement  with  the  experimental  result  might  be  argued.  However  this  may  well 
be  fortuitous  in  that  at  these  highest  frequencies  we  are  approaching  the 
overall  system  noise.  Certainly  the  "too  high"  experimental  results  above 
20  kHz  at  a  =  +30°  are  not  approached  by  the  computed  results  including 
velocity  turbulence. 

The  broad  conclusion  must  be  that  velocity  turbulence  does  appear  to 
have  an  effect  on  the  spectrum  at  the  largest  angles  of  look.  However,  except 
for  -30°,  the  effect  is  within  the  accuracy  of  measurement. 
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Calculation  of  shock  and  separating  streamline  when  p.  <  p 
- - - - -  .  .  -- — .  ,  .  ■  ,  ,  1  ,a 

The  following  equations  determine  the  radius  of  curvature  of  the 

shock  and  r^  of  the  separating  streamline  by  means  of  a  first  order 
expansion  about  the  separation  point: 


r  =  M.  1 
s  Is 


(Y  +  D 


t  ,  (1) 

a 


R  =  x  r 


(2) 


where 


bo  °2  *  b2  '~o 

b2  C1  ‘  bl  C2 


be,  +  b.  c 
o  1 _ 1  o 

b2  C1  "  b1  °2 


(3) 


b  =  sin  <S  cos  6 
o  s  s 


A  +  B  +  c  cot  6 


b .  =  { 1  -  sin  6  cos  6 

s  s 


1 


-A  +  BD  +  cE 


cot  a  }  sin  a 
s  s 


b  =  cos  (a  -  6  ) 
2  s  s 


c  =  F  cot  a 
o  s 


c.  =  (G  +  FD)  cos  a 
1  s 
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c  _  =  II  sin  (o  -  6  ) 

2  s  s 


2 

sec  a 


2  .  2 
2  M,  sin  a 
1 s_ s 

(M .  sin4*  a  -  1 ) 
Is  s 


2  .  2 
2  M„  sin  a 
1  s  s 


1  *  -  M'’s  sin1  o2 


(M,s  -  ’> 


E 


=  1  + 


1 :  L:ur1  M1  s 


(y  +  1)  2 

— * — -  cosec  a 

2  s 


-  1 


2  .  2 

2  Y  sin  a 

_  1  s _ s 

u2  •  2  Y  -  1 

Y  M,  sm  a  — 1 — - — 
'Is  s  2 


and 


Y  M 


1  s 


(M,s  -  ,) 


H  =  y  M, 


M  is  the  Mach  number  in  the  flow  upstream  of  the  shock,  M0 
Is  2s 


Mach  number  in  the  flow  downstream  of  the  shocko  M 


2s 


and  6 


is  the 
are  known  in 


terms  of  M„  and  a  from  the  oblique  shock  relations; 
1  s  s 


M„  and  a  are 
Is  s 


known  in  terms  of  the  pressure  ratio  across  the  shock  at  separation  and  the 
ratio  of  the  ambient  pressure  to  the  stagnation  pressure  ahead  of  the  shock. 
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TABLE  1 


Sample  computer  output 

Microwave  scattering  using  Kolmogoroff  turbulence  factor  with  scattering 


Title 

Geometry 


Transmitter 

Turbulence 

Accuracy 

Jet 


-  ELP:  Variable  A:  Slope/2;  X-band 

-  Receiver  from  exit  plane 
L 

Alpha 
Breakpoint 

-  Wavelength  (Lambda) 


Space  scale  (A) 
Intensity  (i) 


-  G 


Length 

Diameter  at  nozzle 
Maximum  diameter 


Potential  Core  -  Length 

Velocity 
Electron  density 
Collision  frequency 

Angle  of  polarisation  of  zero  with  attenuation 

Carrier  a 1 1  enua t ion 


OoOOO  m 
0.223  m 

0.0873  rad  (5.000  deg) 
1.532  m 

0.0316  m 

0o051  m  (maximum) 

1.000 

12 

3.830  m 
0.153  m 
0.331  m 

1.915  m 

2.290^^+03  m/sec 
2.0001n+l6  m~3 
2.00010+11  sec-1 


-7.573  db 


Bandwidth  30  Hz 


Frequency, 

Hz 

100 

200 

500 

1000 

Bandwidth  500  Hz 

Frequency , 
Hz 

1000 

2000 

5000 

10000 

20000 

50000 

100000 


Noise  power  below  carrier 
per  unit  bandwidth,  db 

-7*1.27 

-75.6*1 

-80.68 

-85.31 


Noise  power  below  carrier 
per  unit  bandwidth,  db 

-85.30 
-91.30 
-100.3 
-108.2 
-118.2 
-  i*+o .  U 


Integrated  power  is  -1*7.28  db  below  incident  carrier 


Apparent  carrier  attenuation  is  -7.571  db 


loss . 
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TABLE  2 


Information  recorded  on  tape 


Channel 

Record 

Module 

Informat  ion 

1 

FM 

Longitudinal  vibration  (RX  aerial) 

2 

DR 

AM/lov  frequency 

3 

FM 

Phase  monitor 

It 

DR 

AM/medium  frequency 

5 

FM 

Signal  AGC 

6 

DR 

AM/high  frequency 

7 

FM 

Reference  AGC 

8 

DR 

PM/low  frequency 

9 

FM 

PSD  output 

10 

DR 

PM/medium  frequency 

1 1 

DR 

Reference  Ail 

12 

DR 

PM/high  frequency 

13 

DR 

Shuttle  pulse  (loop  analysis  purpose 

14 

DR 

Attenuation 

19 

Speech  and  timing  (half  track) 

Alternative  inputs 

avai] able 
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e 
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b) 

dV 

* 

dtt 
L 
r . 


m 


m 


r . 
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transmitter  power 
received  noise  power 
received  signal  power 
transmitting  antenna  gain 
receiving  antenna  gain 

angle  of  sight  at  transmitter  (Fig.  l) 
angle  of  sight  at  receiver  (Fig.  1 ) 
angle  of  scatter  (Fig.  l) 
distance  of  transmitter  from  element  dV 
distance  of  receiver  from  element  dV 

distance  from  receiver  to  transmitter 

-3 

electron  density,  cm 

”3 

fluctuating  part  of  electron  density,  cm 
mean  electron  density,  cm 
refractive  index 

fluctuating  part  of  refractive  index 

mean  refractive  index 

Doppler  frequency  shift ,  Hz 

classical  electron  radius  =  2.8  .  10  cm 

electronic  charge,  esu 

electron  collision  frequency,  sec  ^ 

2tt/X,  cm 

incident  wavelength,  cm 

s  uttering  cross  section/unit  volume,  cm  1 

total  scattering  cross  section/unit  volume,  cm  1 

turbulent  intensity 

turbulent  scales,  cm 

angular  radio  frequency 

spectrum  function 

elementary  volume 

angle  between  direction  of  scatter  and  incident  E  vector 

elementary  solid  angle 

core  length,  cm 

exit  plane  radius,  cm 

radius  of  jet  when  inviscid  boundary  parallel  to  axis,  cm 
distance  from  exit  plane  of  point  where  jet  radius  =  r^,  cm 
core  radius,  cm 
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? 

r 

x 

M 

m 

u 


u. 

1 


P. 

J 


P 

p 

core 

Sc. 


radius  at  which  u  =  u./2,  cm 

1 

radial  co-ordinate,  cm 
axial  co-ordinate,  cm 

centre-line  Mach  number  at  x 

.  .  m  .  -i 

longitudinal  component  of  velocity,  cm  sec 

velocity  at  edge  of  core  if  x  <  L  ,  or  on  the  centre-line  if 

T  -1 

x  >  L,  cm  sec 

centre-line  velocity,  cm  sec  ^ 
exit  plane  pressure 
exit  plane  gas  density 
ambient  pressure 
ambient  gas  density 
gas  density 
core  gas  density 
turbulent  Schmidt  number 
r  -  r. 

l 

(2  ln(2))* 

r i  -  r. 

2  1 


f  mass  fraction  of  original  jet  gas  in  the  mixing  region 

f^  value  of  f  at  the  edge  of  the  core  if  x  ^  L  or  centre-line 

value  of  f  if  x  >  1 

-3 

N  electron  density  m  the  core,  cm 

core  __ 

N  . _  equilibrium  electron  density,  cm 
equil 

v  electron  velocity,  cm  sec 

e 

W  width  of  mixing/shear  layer,  cm 

g  accuracy  parameter 

0  azimuthal  angle  in  cylindrical  co-ordinate  system  within  which  the 

jet  is  defined 
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